Tidal disruption flares from stars on eccentric orbits 

Kimitake Hayasaki'-^'^ ", Nicholas Stone^, and Abraham Loeb^ 

' Department of Astronomy, Kyoto University, Kitashirakawa-Oiwake-cho, Sakyo-ku, Kyoto 606- 
8502, Japan 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA, 02138, USA 
^ Korea Astronomy and Space Science Institute, Daedeokdaero 776, Yuseong, Daejeon 305-348, Ko- 



Abstract. We study tidal disruption and subsequent mass fallback for stars approaching 
supermassive black holes on bound orbits, by performing three dimensional Smoothed 
Particle Hydrodynamics simulations with a pseudo-Newtonian potential. We find that the 
mass fallback rate decays with the expected -5/3 power of time for parabolic orbits, albeit 
with a slight deviation due to the self-gravity of the stellar debris. For eccentric orbits, 
however, there is a critical value of the orbital eccentricity, significantly below which all 
of the stellar debris is bound to the supermassive black hole. All the mass therefore falls 
back to the supermassive black hole in a much shorter time than in the standard, parabolic 
case. The resultant mass fallback rate considerably exceeds the Eddington accretion rate 
and substantially differs from the -5/3 power of time. 

1 Introduction 

There is substantial evidence that galactic nuclei harbor supermassive black holes (SMBHs), the ma- 
jority of which are quiescent and not active galactic nuclei. The tidal disruption of a star by a SMBH, 
and subsequent flaring activity, provides a rare observational diagnosis for the large population of qui- 
escent SMBHs. These powerful flares are expected to have a luminosity at least comparable to the 
Eddington luminosity II1I2I . 

The standard picture of a tidal disruption event (TDE) involves a star at large separation falling 
into a massive black hole on an almost parabolic orbit. After the star is tidally disrupted by the SMBH, 
half the stellar debris becomes gravitationally bound to the SMBH as it loses orbital energy inside the 
tidal radius. The bound debris finally falls back and accretes onto the black hole. Kepler's third law 
implies that the accretion rate decays with the -5/3 power of time IIH3| . 

Observed light curves are in reasonable agreement with this theoretically predicted mass fallback 
rate, although some show deviations |4| and the sample size is sufficiently small to make detailed test- 
ing of theoretical models difficult. Observations suggest that the TDE rate is ~ lO^^yr"' per galaxy ||5l. 
This observed rate is in rough agreement with theoretical rate estimates based on two-body scattering 
at ~ pc scales, which motivates the assumption of nearly parabolic orbits |6|. 

However, recent theoretical studies on rates of tidal separation of binary stars by SMBHs suggest 
that a significant fraction of tidal disruption flares may occur from stars approaching the black hole on 
somewhat eccentric orbits, significantly less parabolic than in the standard picture |7|. Other sources 
of TDEs from stars on more eccentric orbits include binary SMBH systems and recoiling SMBHs |8|. 
These latter two sources are capable of producing TDEs with even lower values of orbital eccentricity 
than in the binary separation scenario, and motivate our work here. In this paper, we explore through 
hydrodynamical simulations how mass fallback rates in TDEs vary between the canonical, parabolic 
case and the underexplored eccentric scenario. 
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2 Numerical method 

We describe here procedures for numerically modeling the tidal disruption of stars on bound orbits. 
The simulations presented below were performed with a three-dimensional (3D) Smoothed Particle 
Hydrodynamics (SPH) code, which is based on a version originally developed by Ref. |9|. We model 
the initial star as a polytropic gas sphere in hydrostatic equilibrium. The tidal disruption process is 
then simulated by setting the star in motion through the gravitational field of an SMBH. 

A star is tidally disrupted when the tidal force of the black hole acting on the star is stronger than 
the star's self-gravity. The radius where these two forces balance is defined as the tidal disruption 
radius 

n = n, (1) 

where Mbh is the black hole mass and m, is the stellar mass. The star-black hole system is put on the x- 
y plane, where both axes are normalized by and the black hole is put at the origin of the system. The 
initial position of the star is given by ro = (ro cos (po, ro sin (po, 0), where ro = 3rt is the radial distance 
from the black hole and ipo shows the angle between x-axis and ro. In our simulations, the black hole 
is represented by a sink particle with the appropriate gravitational mass Mbh- All gas particles that fall 
within a specified accretion radius are accreted by the sink particle. We set the accretion radius of the 
black hole as equal to the Schwarzshild radius rs - 2GMbh/c^, with c being the speed of light. 

In order to treat approximately the relativistic precession of a test particle in the Schwarzschild 
metric, we incorporate into our SPH code the following pseudo-Newtonian potential IfTOl : 



1 - ci rs 

Ci H h C3 — 

l-C2(rs/2r) 2r 



(2) 



where we adopt c I = (-4/3)(2-H V6), C2 = (4 V6-9), and C3 = (-4/3)(2 V6 - 3). Equation © reduces 
to the Newtonian potential when ci - I and C2 = C3 = are adopted. Note that equation (|2]i includes 
no higher-order relativistic effects such as the black hole spin or gravitational wave emission. 

We have performed five simulations of tidal disruption events with different parameters. The com- 
mon parameters through all of simulations are following: m, - IMq, r* = IRq, Mbh - lO^M©, (po = 
-0.47r, and 7 = 5/3. The total number of SPH particles used in each simulation is 10^, and the termi- 
nation time of each simulation is where £3^^ = ^jrl| Gm, ^ 5.1 x lO"^(r»/7?0)-^^^(M0/m»)'^^ yr. 
We also adopt standard SPH artificial viscosity parameters aspH = 1 or jSsph = 2. Table 1 summa- 
rizes each model, where the penetration factor /3 represents the ratio of the tidal disruption radius to 
pericenter distance, rp. 



3 Tidal disruption of stars on bound orbits 

As an approaching star enters into the tidal disruption radius, its fluid elements become dominated by 
the tidal force of the black hole, while their own self-gravity and pressure forces become relatively 
negligible. The tidal force then produces a spread in specific energy of the stellar debris 

GMbh^-. 

Ae ^ . (3) 

The total mass of the stellar debris is defined with the diff'erential mass distribution m{e) = dM(e)/de, 
where M(e) = m(e')de' . When a star is disrupted from a parabolic orbit, m(e) will be centered on 
zero and distributed over -Ae < £ < Ae. 

Since the stellar debris with negative specific energy is bound to the SMBH, it returns to pericenter 
and will eventually accrete onto the black hole. If we define its binding energy, e - -GMbh/^ci 
(the semi-major axis of the stellar debris is a), then the mass fallback rate is given by |[3l dM/dt = 
{dM{£)/de)\d£/dt\ (e < 0), where de/dt = -(l/3)(27rGMBH)^^^r5^^ This is derived from Kepler's 
third law. 
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Table 1. The first column shows each simulated scenario. The second, third, and fourth columns are the penetra- 
tion factor p = r^jrj, the initial orbital eccentricity of star-black hole system e,, and its initial semi-major axis a,, 
respectively. The last column describes the remark for each model. 
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Fig. 1. Differential mass distributions over specific energy of stellar debris and their corresponding mass-fall back 
rates. While the differential mass distribution is shown in panel (a), the mass fallback rate is shown in panel (b). 
In both panels, the dot-dashed line (green), dotted line (blue), dashed line (red), dashed-three-dotted line (black), 
and solid line (black) represent the mass distributions and corresponding mass fallback rates of Model 1 , Model 2, 
Model 3, Model 4 and Model 5, respectively. The energy is measured in units of Ae given by equation l[3j. 



The specific orbital energy of a star on an eccentric orbit is given by 

forb ~ — z = z - e.), (4) 

2a* 2rt 

where a» and e, are the initial semi-major axis of the star-black hole system and its initial orbital 
eccentricity, respectively. This quantity is less than zero because of the finite value of a,, in contrast to 
the standard, parabolic orbit of a star. If eorb is less than /)e, all the stellar debris should be bounded by 
the black hole, even after the tidal disruption. The condition Corb = therefore gives a critical value 
of orbital eccentricity of the star 

^1/3 



below which all the stellar debris should remain gravitationally bound to the black hole. The critical 
eccentricity is evaluated to be Ccrit = 0.98 for Model 4, whereas Ccrit = 0.996 for Model 5. For the ec- 
centric TDEs, the orbital period of the most tightly bound orbit, fmin, and the orbital period of the most 
loosely bound orbit, fmax, are obtained by using Kepler's third law with e - />£+ eo^b and equation (|4]l. 
The duration time of mass fallback for eccentric TDEs with e, < Ccrit is thus predicted to be finite and 

can be written by zff = W-fmin = (tt/ V2)(Q;'/[y6(l-e*)]'^') ([l/2 - (1/^1 - e.)(mJMmy'^l"^ - l). 

Evaluating this gives At 207^3;^' for Model 5, whereas zif — > oo for Model 4 in spite of smaller fmin 
than that of Models 1-3. 

Figure [1] show differential mass distributions and their corresponding mass fallback rates in Mod- 
els 1-5. While the differential mass distribution is shown in panel (a), the mass fallback rate is shown in 
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panel (b). In panel (b), the horizontal solid line denotes the Eddington rate: Mem - (l/'7)(i'Edd/c^) - 
2.2 X 1O"^(77/O.1)"'(Mbh/1O^M0) Mgyr"', where Lndd = 47rGMBH'Wpc/crx is the Eddington luminos- 
ity with nip and ctj denoting the proton mass and Thomson scattering cross section, respectively, and 
77 is the mass-to-energy conversion efficiency, which is set to 0.1 in the following discussion. 

In panel (a), the central peak of Model 1 is attributed to mass congregation, from the self-gravity 
of the stellar debris. The energy spread corresponds to Ae before and after the tidal disruption. The 
corresponding mass fallback rates are proportional to r^^^. The slight deviation from time to the -5/3 
power originates from the convexity around Ae and the central peak rising from Q.2Ae to (see also 
IfTTII ). Simulations of Models 2 and 3 have performed with the pseudo-Newtonian potential given by 
equation (|2]). Model 3 has the same simulation parameters as Model 2 except for 13-5. Since the 
potential is deeper as fi is higher, the re-congregation of the mass due to the self-gravity of the stellar 
debris is prevented. This leads to the mildly-sloped mass distribution, and therefore the peak of the 
mass fallback rate also smooths. 

The mass is not distributed around zero but around -Ae in Model 4, and around —5Ae in Model 5. 
This is because the specific energy of initial stellar orbit is originally negative (see equation]?]). Clearly, 
most of mass in Model 4 is bounded by the negative shift of the center. The resultant energy spread is 
slightly larger than we analytically expected. This suggests that the critical eccentricity is smaller than 
the value in equation (]5]l. In Model 5, all of mass is bounded and falls back to the black hole in a much 
shorter time than that of Models 1-3. As shown in panel (b), the mass fallback rate of Model 5 is four 
orders of magnitude greater than the Eddington rate. 



4 Concluding remarks 

We have performed 3D SPH simulations of tidal disruption processes for stars on bound orbits. Our 
main conclusions are summarized as follows: 

1. There is a critical orbital eccentricity below which all stellar debris falls back to the black hole. 
The simulated critical eccentricity is slightly lower than expected from our analytical prediction. 

2. In an eccentric TDE with orbital eccentricity below the critical eccentricity, all the stellar debris 
falls back to the black hole in a much shorter time than that of the standard TDE. The resultant 
mass fallback rate substantially exceeds the Eddington rate and differs from the -5/3 power of time. 

The full details of this work can be seen in Ref. IHl. 
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